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An  analysis  was  made  <ff  the  optical  recording  and  readout  characteristics  of  the  photodichroic 
KF:LiF  crystal  containing  F^  defects  as  the  activated  medium.  The  results  of  a aeries  of  measuremenU 
indicate,  in  agreement  with  the  analysis,  that  this  crystal  system  has  the  highest  photosensitivity  of  any 
direct  reversible  recording  material  without  a gain  mechanism  (*^1  mJ/ct^.  The  modulation  transfer 
function  for  a 200-|^-thick  activated  layer  was  measured  to  be  40%  at  75  llnes/mm,  and  the  dynamic 
range  exceeded  40  w.  The  crystal  can  be  operated  faUgue  free  at  220  K by  use  of  thermoelectric  ^ 
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OPTICAL  RKCOKDING  CHARACTKRISTICS  OK  KK:l.iK 
PHOTODICHROIC  CRYSTAIJS 


INTRODUCTION 

Backi^round 

A photodk'hruic  material  is  one  which  exhibits  photu-induced  dichroism  (absor|>tion 
of  incident  light  dep<>nils  on  its  |M>lari/ationl  and  birefnngence  (velocity  of  the  light  in  the 
maU'nal  depends  on  its  polarization).  Alkali  halide  crystals  containing  anisotmpic  defects, 
which  can  be  aligned  along  discn>te  crystallographic  directions,  are  well  known  reversible 
photodichroic  materials.  .Much  of  the  previous  work  with  thi's*-  materials  has  dealt  with 
information  storage  and  nondestnctive  retrieval  |1-3|.  The  M and  Ma  defects  | 1.5|  are 
pnme  candidates  for  nondestructive  readout,  but  they  suffer  from  one  or  mor<‘  problems 
such  as  low  photosensitivity,  inconvenient  wavelengths,  or  s<*vere  cooling  n‘<|uir<  ents  to 
avoid  fatigue.  The  P’a  defects  in  KCl:Li('l  and  KCliNaCl  have  also  re«-eiv**d  considerable 
attention  (61.  They  have  the  highest  photosensitivity  of  any  photodichroic,  but  they 
generally  have  a destructive  readout.  F’hotodichroics  have  also  lieen  us«*d  as  holographic 
storage  media  (7|.  Recent  emphasis  on  optical  s|HH'tnim  analysis  of  rasti*r-ri‘cord»*d 
signals  [8,9|  has  caus(>d  renewed  interest  in  photodichniic  materials,  since  the  requinment 
for  a virtual  nondestructive  r«*ad  is  somewhat  n‘lax»*d  In  this  cas<*  the  signal  is  nvord«*d 
by  a modulated  focused  scanning  laser  beam  in  a rasU'r  format  with  frame  timi's  of  a few 
milliseconds.  The  Fourier  transform  can  bt*  obtained  by  uniform  illumination  so  long  as 
the  detector  used  can  analyze  the  sp«*ctrum  before  significant  decay  of  the  recordi-d 
raster.  It  may  also  be  necessary  to  recover  the  rc‘cordt>d  raster  from  the  input  transducer 
after  analysis  by  re.scanning,  and  this  must  be  accomplisht'd  before  the  raster  is  erased. 
Alternative  input  transducers  such  as  the  General  Klecthc  cohenmt  light  valve  have  no 
storage  facility  and  cannot  accomplish  this.  In  addition  to  a photodichroic  mat4>rial  lieing 
used  as  an  input  transducer  to  an  optical  spectrum  analyz<‘r,  it  has  an  important  poUmlial 
application  in  the  analysis  of  these  signals  in  the  transform  plane  [lO]  In  the  spectral 
analysis  of  the  recorded  raster  the  presence  of  persistant  strong  narrow-band  sources  makes 
it  difficult  to  detect  adjacent  weak  sources  and  particularly  transient  sources.  The  photo- 
dichroic  material  can  operate  as  a change  detector  in  this  case  or  attenuate  lung-U>rm 
strong  background  signals  while  transmitting  weak  adjacent  signals.  (This  application, 
known  as  signal  excision,  is  described  in  Appendix  A.)  A schematic  of  both  uses  of  the 
photodichroic  is  shown  in  Fig.  1. 

This  report  describes  developmental  work  on  a relatively  new  photodichroic  crystal: 

Fa  defects  in  KF:LiF.  Since  the  crystal  is  hydroscopic,  special  handling  is  required,  but 
it  has  many  desirable  features  not  shared  by  other  photodichroics.  Its  photosensitivity  is 
among  the  highest  of  all  photodichroics,  it  can  be  operated  fatigue  free  not  far  below  room 
temperature  using  simple  thermoelectric  coolers  in  a dry  nitrogen  atmosphere,  and  its 
wavelength  range  is  ideally  matched  to  the  strong  argon  laser  lines. 


Manuscript  aubmitted  Auipiit  26,  1977. 
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Fig.  1 — Raster  scanned  (iptical  spectral  analyzer  using  a photodichroic  as  an  input  transducer  and  as  an 
adaptive  spatial  filter  in  the  transform  plane.  The  raster  scanner  is  shown  off  axis,  but  it  could  be  placed  on 
axis  by  use  of  a beamsplitter. 


Physics  <»f  the  Fa  Defect 

When  a normally  transparent  alkali-halide  crystal  is  exposed  to  ionizing  radiation  or 
heatsKl  in  an  alkali  atmosphere,  halide  vacancies  are  formed  which  after  trapping  an  elec- 
tron can  absorb  visible  li^t.  The  crystal  is  thus  colored,  and  the  defects  formed  are 
called  F centers.  The  absorption  of  F centers  is  isotropic  and  occurs  in  a single  Gaussian 
band  with  no  useful  properties  as  a recording  medium.  However,  if  the  crystal  contains 
alkali  impurities  of  smaller  ionic  radius  than  the  host  (such  as  the  KF;LiF  mixture)  and 
the  F centers  are  optically  excited  at  the  appropriate  temperature,  the  defects  diffuse  and 
become  trapped  at  sites  adjacent  to  the  Li*  sites.  The  new  defects  formed  are  known  as 
F^  centers,  and  their  optical  properties  are  radically  altered  (11,12].  Fa  centers  (Fig.  2) 
are  anisotropic,  with  a triad  of  optical  dipole  moments,  one  being  parallel  to  the  vacancy- 
impurity  axis  and  two  being  perpendicular  to  this  axis.  The  usual  F-center  singlet  absorp- 
tion splits  and  becomes  a doublet  corresponding  to  the  triad  of  optical  moments  along 
the  crystal  < 100  > axis.  The  moment  parallel  to  the  impurity-vacancy  axis  gives  rise  to 
a longer  wavelength  absorption  (Fa,  ),  and  the  two  perpendicular  moments  are  degenerate 
and  give  rise  to  a single  shorter  wavelength  absorption  (Faj).  All  three  moments  have 
equal  oscillator  strengths. 

Thus  the  Fa  absorption  consists  of  a double  band  whose  components  partially  over- 
lap. In  addition,  when  these  Fa  defects  are  optically  excited,  they  can  rotate  about  the 
Li*  defect  by  ionic  motion  of  the  halide  ions  to  a new  perpendicular  orientation  with  a 
quantum  efficiency  of  0.5.  As  long  as  the  Fa  centers  have  their  impurity-vacancy  axes 
randomly  oriented  (1/3  along  each  of  the  three  < 100  > directions),  the  crystal  absorbs 
light  isotropically.  However  if  the  Fa  centers  are  excited  with  light  in  one  band  exclusively 
and  polarized  along  the  crystal  axis,  the  defects  rotate  until  their  optical  moment  is 


NRL  RKPORTMIfiH 


0 


© 


© 


© 


(OO.) 


(OlO] 


(.ool 


Fig.  2 — The  defect  in  a KF  LiF  lattice.  Thia 
conaiat  of  a F"  vacancy  which  haa  captured  an  elec- 
tron, with  a Li'*^  nearest  neighbor.  The  three  optical 
dipole  momenta  are  shown 


perpendicular  to  the  polarization.  Thus  the  centers  can  be  preferentially  ali|pti>d.  The 
medium  acts  as  a temporary  polarizer  to  incident  light  in  either  of  the  absorption  bands. 
The  polarizing  property  is  lost  by  randomizing  the  defect  orientations,  either  by  exposing 
to  light  whose  polarization  bisects  the  crystal  axes  or  using  a wavelength  which  equally 
excites  both  absorption  bands.  Thus  the  medium  behaves  as  a reversible  polarizer  whose 
polarization  properties  can  be  spatially  and  temporally  controlled  by  the  polarization  of 
incident  light.  T^ese  properties  are  best  used  when  the  crystal  is  placed  between  cross<*d 
polarizers  whose  easy  directions  are  45°  to  the  crystal  axis  (Fig.  3).  In  those  art'as  of  the 
crystal  where  the  defects  are  randomly  oriented,  the  crossed  polarizers  extinguish  the 
readout  light.  For  nonrandom  orientation  the  polarization  direction  of  the  readout  light 
is  rotated  in  penetrating  the  crystal,  and  the  analyzer  transmits  light  from  that  area.  The 
transmission  is  determined  by  the  relative  absorption  of  light  along  the  crystal  axis,  that 
is,  by  the  dichroism  induced  by  nonrandom  orientation. 


THEORETICAL  ANALYSIS 


Signal  Recording  Characteristics 


Signals  are  recorded  on  a photodichroic  crystal  by  optically  pumping  anisotropic 
defects  from  a random  distribution  to  an  aligned  distribution  as  described  in  Appendix  A. 
The  efficiency  for  inducing  alignment  is  proportional  to  the  difference  in  the  fraction  of 
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Kik  3 — Photodirhruir  rrcordinn  and  rradinR  uperatiuna.  The  light  dbtributiun  in  Ihr 
write  tx’am  modulatea  the  pulari/atiun  prupertiei  of  thr  pbotodichroir  crystal  The 
uniriirm  readout  beam  has  its  polariastinn  rotated  in  thr  same  pattern  as  the  srritr 
beam  Only  the  rotated  parts  of  the  readout  beam  are  thus  transmitted  by  the  analyzer. 
The  output  light  distribution  duplicates  thr  write-beam  distribution  Thr  off  axis 
recording  lieam  could  be  placed  on  axis  by  use  of  a beamsplitter. 


intid»>nt  light  ({wjlarized  along  the  crystal  axis)  absorbed  by  the  components  Fa,  and  Faj. 
If  the  wavelength  is  such  that  both  bands  are  equally  excited,  then  no  alignment  occurs. 
The  wavelength  dependence  of  the  recording  sensitivity  for  a KF:I,iF  crystal  at  240  K is 
shown  in  Fig.  4. 

The  dynamic  characteristics  of  the  recording  mtKlium  can  be  calculated  for  exposure 
U)  light  of  arbitrary  wavelength  and  polarization.  This  involves  three  coupled  nonlinear 
differential  equations  which  are  quite  complicated.  However,  for  the  simple  and  important 
case  of  excitation  in  one  band  exclusively  and  with  polarization  along  the  crystal  axis,  the 
differential  equations  simplify  and  are  derived  and  solved  in  Appendix  B.  The  solutions 
for  the  optical  densities  measured  along  the  crystal  axis  are  given  by  Fqs.  B12  as 

K,(E)«logI(10^“  - De-Pf"'  + 1 

and 

KiiE)-  I /fo  - I KiiE} 

where  Kq  is  the  initial  optical  density  measured  along  either  crystal  axis  (assuming  initisd 
randomization),  E is  the  exposure  in  mJ/cm^,  andp  is  a factor  in  cm^/mJ  depending  on 
temp<‘rature  as  given  in  Appendix  B.  Thus,  whe.i  a .signal  is  recorded  onto  a photodichroic, 
it  is  stored  as  a difference  in  absorption  along  the  crystal  axis;  Ky  - K^.  These  calculated 
recording  characteristics  are  plotted  in  Fig.  5. 
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nm  for  ■ KF  LiP  crytui  at  240  K and  an  initial 
optical  drnaity  of  0.6.  The  aicnal  ii  recorded  aa  a 
difference  in  optical  denaity  alone  the  crytlal  axia 
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SiKnul  Reo»v»'r>-  Charac-tensUi^ 

rh«‘  siKiial  IS  rit'Dvoriti  by  unifonniy  illuminatinK  tht*  crystal  with  littht  (xilanzi'd  at 
45  to  the  crystal  axis.  ITii*  li|{ht-int(>nsity  dtstnhutiun  transmitted  through  an  analyzi'r  on 
the  output  side  of  the  crystal  duplicaU's  the  oriKinal  r<‘<'ordin|{  lit;ht  distnhution,  as 
s<-hematically  shown  in  Kik  3.  The  transmi.ssion  aris«‘s  from  the  s«‘UH'tive  attimuation  and 
rotation  of  the  incident  t'  vtH-tor  comfxments  aJon^  the  crystal  axis.  The  transmission 
e<|uation  IS  denved  in  .Apia-iulix  ('  as  Kq.  ('H  and  is 

r-  ^ ♦ 10  **''**‘*''2sin2  &/2. 

where  & is  the  phasi-  anitle  induced  hy  the  hiiefnnitence  given  by  Kqs.  (’2  and  (?3. 

This  (‘<]uation  stat4*s  that  the  absorptive  com|M>nent  of  the  transmission  cannot  ex- 
ce«*d  25'/f  hut  the  phas«-  component  can  approach  100'?  for  sufficiently  low  absorption 
and  large  phas<‘  angle.  Of  coursi*  iMith  of  these  b-rms  cannot  be  a maximum  .simultani'ously. 
Kigun*  6 shows  the  wavelength  de|H‘ndence  of  the  transmission  function  (solid  curve  in  the 
lower  part  of  the  figure)  for  a saturation  write  indicated  by  the  dichroic  ahsorfition  in  the 
iq>p«“r  jiart  of  the  figure.  The  absorjitive  and  phase  contributions  to  the  transmission  func- 
tion are  shown  by  dottinl  and  dash«>d  curves  res|H*ctively.  The  s<*ven*  atUmuation  in  trans- 
mission lielow  450  nm  is  due  to  the  higher  average  absor|>tion  in  this  wavelength  range. 

The  transmission  at  48H  nm  is  purely  due  to  |>hasi>  effects,  but  at  514  nm  there  is  a nearly 
iKjual  contribution  from  the  absorptive  term.  The  dichroism  shown  in  Fig.  6 was  chos«*n 
atypically  large  to  illustrate  the  {xiU'iitially  large  transmission  available  at  488  nm.  In 
Fig.  7 a s«*ri»*s  of  tran.smis.sion  curvi*s  are  shown  coire.sponding  to  sucessive  changes  in 
dichroi.sm  as  indicated  by  the  scale  changes.  The  phas(‘  contribution  grows  faster  than  the 
absoriitive  for  larg«-r  dichroism.  Thesi'  transmission  characteristics  are  for  a crystal 
U-mjx-rature  of  240  K.  Since  the  dichroii-  bands  shift,  become  more  narrow,  and  increast* 
in  height  as  the  U-mfH-rature  is  lowen-d,  the  transmi.ssion  curves  will  change.  The  dichroic 
absorption  at  77  K is  shown  in  the  up|><>r  part  of  Fig.  8 along  with  the  com*s|H)nding 
transmission  cari.’e  b«*low.  The  two  dichroic  bands  an-  now  almost  fully  s»*parat4*d,  and 
the  average  optical  density  at  488  nm  is  quite  small.  This  accounts  for  the  large  incrf>as** 
in  transmission  at  this  wavelength  over  that  observ«-d  in  the  higher  U‘m|H*ratun’  case  in 
P’ig.  6.  Figure  9 .shows  the  effwts  of  varying  the  dichroism  scale  on  transmi.ssion  at  77  K. 
The  transmission  at  488  can  become  quite  large  for  the  higher  densities.  The  transmi.ssion 
is  still  significant  out  to  nearly  600  nm.  This  high  transmission  corresponding  to  large 
dichroism  is  obtairM'd  at  the  expens«>  of  n>solution,  as  will  b«*  explained  later. 


Dynamic  (Characteristics 


The  complete  dynamic  behavior  for  writing  at  514  nm  can  be  obtained  by  substitut- 
ing PCqs.  B12  into  Kq.  C8.  With  the  appropriate  parameters  for  temperature  and  initial 
optical  density  this  gives  the  transmission  at  any  wavelength  as  a function  of  the  laser 
energy  exfrosure.  The  results  of  these  calculations  are  shown  in  the  series  of  log-log  plots 
of  Figs.  10a  through  lOd.  The  changes  in  the  Kq  values  labeled  on  the  curves  from  the 
240-K  to  the  77-K  plots  are  due  to  the  narrowing  and  increased  height  of  the  bands.  This 
accounts  in  part  for  the  enhanced  sensitivity  shown  at  77  K.  F'or  low  exposures  all  of 
these  curves  are  approximately  linear  in  amplitude  transmittance.  However  there  are 
.significant  dynamic  differences  Ix'tween  the  514  nm  and  488-nm  readouts  for  changes  in 
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Pig  6 — Rpidiiut  tranimiacion  rhainrtrrictioii  at  240  K for  tha  photo 
inducpd  dirhroiam  ihown  at  tha  tup  of  tha  figura  fora  514  nm  wnla 
to  laturation  Tha  daahad  curva  and  tha  duttad  rurva  giva  tha  phaaa 
and  tha  absorptiva  mmponanta  of  tha  tranammion  raapactivaly. 

the  parameter  Kq.  'rhis  is  most  evident  in  the  514-nm  readouts  at  77  K (Kir.  10c).  The 
slopes  chanRe  as  the  value  of  Kq  is  chanRed,  in  contrast  to  the  488-nm  curves  (FiR.  lOd). 
The  origin  of  this  behavior  lies  in  the  difference  in  dynamic  behavior  between  the  phase 
and  absorptive  components  of  the  transmission.  Since  the  514-nm  readout  is  composed 
of  a mixture  of  these  effects,  its  dynamic  characteristics  are  more  complicated. 

The  difference  in  dynamic  behavior  between  the  514-nm  and  488-nm  readouts  is 
illustrated  by  plotting  d \/TldE  vs  exposure.  For  maximum  linearity  these  functions 
should  be  nearly  constant  with  exposure.  The  results  of  the  calculations  are  shown  in 
Figs.  11a  through  lid.  These  graphs  illustrate  more  clearly  how  the  linearity  and  sensi- 
tivity change  with  changes  in  the  initial  optical  density  Kq.  The  linearity  changes  with 
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EXPOSURE  (mJ/cm*) 

(b)  488-nm  read  at  240  K 

Pi*.  10  — Readout  transminion  vs  input  exposure  calculated  from 
Eqs.  B12  and  C8  for  814-nm  write  and  the  indicated  crystal  tempera- 
ture and  read  wavelength.  The  upper  straight  line  indicates  a slope  of 
0.5  in  each  log-log  plot,  which  corresponds  to  an  ideal  square-law 
response  (Continued). 
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Fig.  10  — Readout  tranimiuion  vs  input  exposure  calculated  from 
Elqa.  B12  and  C8  for  5I4-nm  write  and  the  indicated  crystal  tempera- 
ture and  read  wavelength.  The  upper  straight  line  indicates  a slope  of 
0.5  in  each  log-log  plot,  which  corresponds  to  an  ideal  square-law 
response  (Continued). 
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Fig.  1 1 — Derivative  of  the  amplitude  tranamittance  va  exposure,  ralrulated  from  Eqs.  B1 2 and 
C8  fur  514-nm  write  at  thf  indicated  crystal  temperature  and  read  wavelength  (Continued) 
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Kq  in  all  of  the  plots,  but  for  the  514-nm  read,  certain  values  of  Kq  can  b«‘  (■hus(‘n  to 
achieve  a nearly  flat  response  at  only  a small  sacrifice  in  sensitivity.  This  is  particularly 
true  at  the  lower  temp«'ralure.  Thus  the  plots  of  P'ig.  1 Id  pn>dicts  that  the  curve  for  an 
approximaU'  optical  density  of  1 gives  the  greatest  linearity.  Since  the  origin  of  this  effect 
lies  in  the  admixture  of  the  absorptive  term  in  the  transmission,  similar  results  could  be 
achieved  by  small  shifts  in  wavelength  of  n-adout. 


EXPKRI.MENTAL  RESULTS 
Sen.sitivity 

Measun’ments  of  nvording  sensitivity  are  in  close  agreement  with  the  theoretical 
analysis.  The  response  generally  saturates  at  exposures  of  10  to  20  mJ/cm^  for  moderate 
optical  densities  of  0.4  to  0.6.  The  minimum  exposure  for  a deUvtable  signal  is  limited 
by  the  system  noise,  due  to  optical  im|>erfections  in  the  crystal  and  components  as  well 
as  electronic  noise.  Signals  have  been  ri>covered  for  exposures  as  low  as  0.14  mJ'cm^ 
with  a signal-to-noise  ratio  of  approximately  2.  Since  much  of  the  optical  noise  is  due  to 
crystal  surface  imperfections,  residual  strain,  and  the  limited  extinction  coefficient  of  the 
polarizers,  it  appears  primarily  as  low-frequency  noise  in  the  optical  transform  plane. 

Thus,  when  the  Fourier  transform  of  the  stored  .signal  is  observed,  much  of  the  noise  is 
removed  from  the  desired  signal.  Under  these  conditions  an  expo.sure  of  60  pJ/cm^  gives 
a detectable  signal.  With  sufficient  optical  quality  this  value  should  also  apply  to  signals 
recovered  from  the  image  plane. 

The  quantum  efficiency  for  reorientation  is  independent  of  temperature;  but  as  the 
temperature  is  lowered,  the  dichroic  absorption  bands  shift,  narrow,  and  increase  in  height. 
The  photosensitivity  increases  at  514  and  488  nm  by  approximately  a factor  of  2,  as  was 
predicted  by  the  theoretical  analysis.  Unfortunately  the  only  available  method  of 
achieving  these  low  temperatures  (below  100  K)  is  the  use  of  cryogenic  fluids.  The  boil- 
ing fluid  introduces  vibrations  which  prevent  the  measurement  of  other  critical  parameters 
such  as  the  modulation  transfer  function  and  the  dynamic  range.  The  use  of  a closed- 
cycle  refrigeration  system  would  avoid  these  difficulties.  Even  though  the  quantum 
efficiency  is  temperature  independent,  at  sufficiently  high  temperature  the  relevant  defects 
began  to  thermally  reorient  and  near  room  tem[ierature  this  limits  the  time  to  which 
signals  can  be  stored  to  a day  or  less.  A more  serious  effect  of  high-temperature  operation 
lies  in  the  thermal  diffusion,  which  causes  the  relevant  defects  to  aggregate  and  thus  lose 
their  dichroism.  This  loss  of  dichroism  due  to  diffusion  occurs  at  room  temperature  after 
approximately  10'’  write-and-read  cycles,  but  at  240  K this  loss  doer  not  occur  until  10” 
cycles  or  better.  There  is  some  fluctuation  in  this  parameter,  depending  on  crystal  purity 
and  defect  concentration,  but  generally  a temperature  of  220  K or  lower  assures  operation 
that  is  virtually  fatigue  free. 


Modulation  Transfer  Function 

Since  the  defects  involved  in  recording  are  atomic  dimensions,  the  medium  is  grain- 
less and  the  resolution  is  inherently  diffraction  limited.  However,  when  used  with  low-f- 
number  optics,  the  thickness  of  the  activated  volume  must  not  greatly  exceed  the  depth 
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CRYSTAL  THICKNESS 


Fig.  12  — The  recording  of  two  focused  spots  of  light 
separated  by  the  spot  diameter  (Raleigh  criterion  for 
resolution)  inside  a crystal  when  the  crystal  thickness 
exceeds  the  depth  of  focus.  Even  though  the  beams 
overlap  inside  the  recording  medium,  the  spots  can 
still  be  resolved  on  readout  due  to  the  higher  writing 
intensity  inside  the  depth  of  focus,  although  with  a 
lower  modulation  transfer  function.  The  diagram  is 
highly  exaggerated  for  illustration. 


of  focus  of  the  optics.  For  an  f/8  system  with  a spot  size  of  5 pm  (100  lines/mm)  the 
depth  of  focus  is  approximately 

d = 2(f  number)(5  pm)  = 80  pm.  (1) 

If  the  crystal  thickness  exceeds  this,  the  modulation  transfer  function  (MTF)  will  de- 
crease due  to  the  overlap  of  the  light  cone  away  from  focus,  as  shown  in  Fig.  12.  The 
MTF  can  still  be  significant  even  for  large  variations  from  Eq.  1 due  to  the  higher  record- 
ing intensity  at  the  waist,  where  there  is  no  overlap. 

A typical  MTF  curve  is  shown  in  Fig.  13  for  a crystal  with  an  active  thickness  of 
200  pm.  A depth  of  focus  of  200  pm  corresponds  to  a spot  size  of  12  pm,  or  a resolu- 
tion of  40  lines/mm.  The  MTF  should  begin  to  drop  at  this  frequency,  and  Fig.  13 
confirms  this.  The  drop  is  rather  rapid  between  40  and  70  lines/mm,  but  there  is  a long 
slowly  decreasing  tail  out  to  beyond  120  lines/mm,  which  corresponds  to  the  high- 
intensity  recording  at  the  beam  waist. 

These  measurements  were  made  by  several  techniques.  The  simplest  method 
involved  recording  the  image  of  a bar  chart  onto  the  crystal  and  recovering  the  stored 
image  onto  photographic  film  after  40X  magnification.  The  MTF  was  obtained  from  the 
film  with  a microdensitometer.  A second  technique  involved  recording  an  8-to-12-Iine/mm 
Ronchi  grating  onto  the  crystal  and  observing  the  Fourier  transform  of  the  stored  grating. 
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Ki((.  13  — Modulation  transfer  function  of  a KF  LiF  photodichroic 
crystal  with  a rccordinn  thickness  of  200  pm 


By  comparing  Iht*  relative  magnitude  of  the  diffraction  orders  with  that  of  the  original 
grating,  the  MTF  is  obtained.  A third  technique  involved  exposing  the  crystal  to  a 
scanning  modulated  laser  beam  and  recovering  the  stored  signal  with  uniform  illumination 
onto  photographic  film.  The  MTF  was  obtained  from  the  film  with  a microdensitometer. 


Linear  Dynamic  Range 

The  dynamic  range  is  limited  for  large-signal  recording  by  nonlinearity  and  for  small- 
signal  recording  by  the  noise  generated  by  surface  imperfections  and  bulk  strain.  Since 
the  material  fabrication  is  still  in  an  early  stage  of  development,  the  optical  quality  of  the 
.surfaces  and  the  strain  were  such  that  noise  prevented  the  detection  of  nonlinearity  for 
the  lowest  exposures  (below  2 mJ/cm^).  Hence  the  dynamic-range  data  presented  are  a 
measure  of  the  lower  limit  for  low  exposures  and  are  more  accurately  a measure  of 
linearity  for  the  higher  exposures. 

A sinusoidal  grating  of  30  lines/mm  was  recorded  onto  a crystal  by  interfering  two 
collimated  laser  beams  at  514  nm.  The  Fourier  transform  of  the  recorded  grating  was 
detected  by  a vidicon  with  a line-scan  output.  Nonlinear  recording  is  indicated  by  the 
presence  of  a second  order  in  the  transform  plane.  By  measuring  the  ratio  of  the  intensity 
of  the  first  order  to  that  of  the  second  order,  the  linear  dynamic  range  in  the  transform 
plane  is  obtained.  Two  sequences  of  measurements  were  performed.  One  involved  using 
nearly  full  modulation  ( %80%)  and  varying  the  basis  level,  and  the  other  involved  varying 
the  bias  but  keeping  the  modulation  fixed  at  0.4  mJ/cm^.  The  results  of  the  first  meas- 
urements are  shown  in  Fig.  14a.  The  second  order  was  below  the  noise  in  the  data  for 
exposures  below  2 mJ/cm^;  hence  a 40-dB  linear  dynamic  range  must  be  considered  a 
lower  limit.  With  improved  surface  conditions  and  strain-free  crystals  the  dynamic  range 
should  exceed  40  dB.  The  linearity  decreases  monotonically  for  increasing  exposure  with 
80%  modulation.  Also  shown  is  the  signal  amplitude  transfer  function  (from  Eqs.  B12 
and  C8),  indicating  the  degree  of  modulation.  In  the  second  measurement  the  bias  level 
was  varied  but  the  modulation  was  fixed  at  0.4  mJ/cm^.  These  results  are  shown  in  Fig. 
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Fig.  14b  — The  measured  linear  dynamic  range  for  a KK:LiF 
photodichroic  as  a function  of  recording  exposure  using  a fixed 
modulation  of  0.4  mJ/cm^.  Also  shown  is  the  amplitude 
transmittance  from  Fig.  10a,  indicating  the  modulation  and 
bias. 


sources  for  KF:LiF  are  NRL’s  Electronic  Materials  Technology  Branch,  Cornell  University’s 
crystal-growth  facility,  and  The  University  of  Utah’s  crystal-growth  facility.  The  latter 
source  has  produced  the  most-strain-free  crystals  of  the  highest  optical  quality  to  date 
and  is  the  major  supplier.  The  usual  cross  section  of  these  single-crystal  boules  is  approxi- 
mately 6 cm2,  but  they  can  be  grown  much  larger. 

Upon  receipt  the  single<rystal  boules  are  carefully  cleaved  (to  avoid  .straining)  into 
slabs  approximately  3 mm  thick.  Since  the  KF;LiF  crystals  are  moderately  radiation 
hard,  it  is  difHcult  to  introduce  a high  concentration  of  F^  defects,  and  special  treatment 
is  required:  an  anneal  of  the  3-mm  slabs  at  700°C  in  a mixed  hydrogen  and  K* 
atmosphere.  This  treatment  causes  H~  ions  to  diffuse  into  the  lattice  and  substitute  for 
the  F'  ions.  The  concentration  of  H"  ions  depends  on  the  hydrogen  pressure  (normally 
5 atmospheres),  and  the  penetration  depth  depends  on  the  temperature  and  annealing 
time.  A 4-hour  anneal  at  700°C  gives  a penetration  of  approximately  200  pm.  In  this 
surface  layer  the  crystal  composition  is  KF:LiF:KH.  After  the  anneal,  the  crystals  must 
be  cooled  under  precise  control  (10° C/mm).  Thermal  strain  is  introduced  if  the  cooling 
is  too  fast,  and  phase  separation  of  the  KF:LiF  mixture  occurs  if  the  cooling  is  too  slow. 
After  the  3 mm  slabs  are  cooled,  they  are  cleaved  so  that  only  one  surface  is  composed 
of  KF:LiP:KH.  Then  the  crystals  are  polished  and  are  sealed  between  glass  flats  using  an 
optical  grade  silicone-rubber  compound  to  protect  the  surfaces  from  moisture.  This  latter 
technique  has  been  used  for  convenience  in  the  early  developmental  stage  but  introduces 
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multiple  n-n«H-tions  due  to  index -matehinii  problems.  A sup*‘nor  Un-hnique  would  b<-  to 
sputU'r  Sj02  ur  evaporaU*  MKF2  films  of  0.5  to  0.1  pm  depth.  Preliminary  t4*8ts  indicate 
the  feasibility  of  achievinK  layers  of  this  thickness,  provided  the  crystal  surfacini  an*  of 
sufficiently  (lood  quality. 

Once  the  crystals  an*  polished  and  sealed,  they  an*  exposed  to  ioni/inx  radiation. 

The  substitutional  H"  ions  are  easily  driven  out  of  their  lattice  siU*s,  and  they  combine, 
forming  stable  H2  molecules.  Tlu‘se  molecul«*s  are  outgas.s<*d  from  the  crystal,  leaving 
b«*hind  negative-ion  vacancies  which  trap  electrons  ( F centers).  The  final  sU*p  involvi*s 
exposing  the  crystals  to  blue  light  near  room  t«*mp«*ratun*.  This  causes  the  F cenU*rs  to 
diffuse,  until  they  become  trap|H*d  at  Li*  sit»*s,  thus  forming  the  relevant  Fy^  centers. 

The  crystals  an*  now  n*ady  to  be  used  aftt*r  cooling  to  at  least  -30’C  to  avoid  fatigue 
problems.  They  are  currently  being  cooled  with  a two-stage  thermoele<  tric  cooler  aft«*r 
l)eing  sealed  in  a dry  nitrogen  atmosphere.  The  requin*d  cooling  is  close  to  the  limit  for 
this  tyj)e  of  cooler.  However  optimum  performance  in  t4*rm.s  of  sen.sitivity  and  dynamic- 
range  will  occur  near  100  K.  Hence  a clos«*d-cycle  n*frigeration  system  would  lie 
requin*d.  Sensitivity  should  improve  by  a factor  of  2,  and  the  dynamic  range  should 
improve  to  bc*tter  than  50  dB.  The  analysis  of  p<*rformance  at  low  tem|M*ratun‘  is  under 
way. 


An  im))ortant  an*a  of  mat«*rial  im|)rovement  involves  increasing  the  conc-entration  of 
defiH  ts.  This  would  allow  thinner  layers  of  activated  defects  and  thus  allow  higher 
resolution.  Currently  the  maximum  concentration  of  stable  defects  achieved  is  approxi- 
mately 10  to  15  parts  per  million.  It  should  he  possible  to  increase*  this  by  several  factors 
h«*fore  exc«*eding  the  solubility  limit  of  LiF  in  KF. 


CONCLUSIONS 

The  characteristics  of  the  KF:LiF  photodichroic  crystal  make  it  a leading  candidate 
for  use  as  an  input  transducer  or  adaptive  spatied  filter  in  an  optical  s|>ectrum  analyzer. 

It  is  less  sensitive  than  the  electrooptic  photoconductors  but  has  much  betU'r  r«‘solution. 
Even  so  the  sensitivity  is  adequate  for  real-time  recording  using  currently  available  argon 
lasers.  Since  the  material  is  still  in  an  early  stage  of  development,  improvements  an*  to 
be  expected  in  crystal  quality  which  will  extend  the  dynamic  range  lx*yond  the  present 
40  dB,  which  is  adequate  for  many  applications.  Since  the  material  nvords  with  maxi- 
mum quantum  efficiency,  no  improvement  in  the  inherent  .sensitivity  can  be  anticipated. 
However,  by  use  of  low-temperature  housing,  the  relevant  bands  can  be  made  to  shift  for 
maximum  use  of  the  514-nm  argon  line,  thus  improving  the  photosensitivity  by  at  least 
a factor  of  2.  A second  area  of  improvement  lies  in  the  increase  in  concentration  of  the 
Fa  defects,  which  will  allow  thinner  coloration  layers.  This  will  improve  resolution,  and 
a modulation  transfer  function  of  50%  at  150  lines/mm  appears  feasible. 
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Appendix  A 

SIGNAL  EXCISION 


ITio  iharacteristics  of  th<*  KF;LiF  photcxlirhraic  when  us<h1  as  a signal  exrisor  in 
the  frequency  plane  of  an  optical  spixtrum  analyzer  will  lx*  described  in  this  ap|>endix. 
Such  an  arrangement  would  apfx*ar  as  shown  in  Fig.  Al.  In  an  optical  spwtrum  analyzer 
the  frequency  .s|MH  trum  of  an  input  U*mporal  signal  ap|M*ars  as  spots  of  light  in  the  trans- 
form plane.  In  the  analysis  of  wide-hand  signals  the  pres«*nce  of  persisU*nt  strong  narrow- 
l>and  soun*es  makes  it  difficult  to  deU'ct  adjacent  weak  oni*s.  The  photodichroic  will  act 
as  a change  deUx  tor  in  this  situation  and  will  atU*nuaU*  long-term  strong  sources  while 
transmitting  weak  ont's. 

Although  there  an*  many  possible  configurations,  the  arrangement  shown  in  Fig.  Al 
is  typical.  Here  the  crystal  is  placed  in  th«*  transform  plane,  and  the  dark  sp<»ts  along  the 
X axis  on  the  crystal  illustrate  the  presence  of  five  signals  of  various  intensities.  These 
light  spots  are  polarized  at  45°  to  the  crystal  axis.  The  indicaUxl  control  light  would  be 
any  broadband  blue  source  which  is  vertically  polarized.  A typical  secjuence  might  be  as 
foi'ows:  The  control  light  source  uniformly  polarizes  the  crystal  with  the  signal  from  the 
input  off.  At  f = 0 the  input  signal  is  allowed  to  fall  on  the  crystal  as  shown  in  F’ig.  Al. 
At  this  instant  all  signals  are  transmitted  with  no  change  in  their  relative  intensity.  How- 
ever, since  the  signal  polarization  bisects  the  crystal  axes,  there  is  a randomizing  effect 
and  thus  a loss  of  dichroism  in  these  exp<«>«‘d  areas.  The  transmission  decay  is  purely  a 
function  of  the  cumulative  energy  of  exposure,  so  the  area  expos«*d  to  the  central  intense 
light  decays  much  more  rapidly  than  the  weakly  exposed  areas.  After  a certain  time  t, 
which  de|)ends  solely  on  the  peak  intensity,  the  strong  signal  is  suppressed  much  more 
si*verely  than  the  weaker  ones,  as  shown  in  the  reimagixl  output  plane. 

'Fhe  precise  dcx*ay  characteristics  depend  on  the  initial  optical  density,  the  initial 
degree  of  alignment,  the  temp<*rature,  and  the  wavelength,  but  a typical  set  of  decay 
curves  is  shown  in  Fig.  A2.  If  the  crystal  is  exposed  to  uniform  light  with  the  indicated 
inU-nsities,  the  transmission  of  light  from  that  area  will  decay  as  shown.  In  the  insert  a 
Gaussian  beam  profile  is  shown.  If  the  crystal  is  exposed  to  such  a beam  profile,  the 
pi*ak  will  decay  as  curve  A and  the  annular  ring  at  B will  decay  as  curve  B.,  etc.  out  to 
the  low-inU*nsity  wings,  'fbus  after  a certain  exposure  such  a Gaussian  beam  would 
ap(M>ar  as  a weak  annular  ring  of  light. 

'Hie  magnitude  of  the  final  decrease  in  transmission  is  ultimately  limited  by  the 
extinction  coefficient  of  the  crossed  polarizers,  which  can  be  from  lO*’’  to  10'®.  After 
a few  milliseconds  of  exposure  (for  a peak  intensity  of  20W/cm2),  the  signal  source  is  cut 
off,  the  control  light  source  realigns  the  crystal,  and  the  system  is  ready  to  be  rec..cled. 
Close  to  a 1009?.  duty  cycle  is  possible,  if  both  signal  and  control  light  .sources  are  of 
sufficient  intensity. 
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F’ig  AI  — Illuxtration  of  the  photodichroic  crystal  used  as  a signal  excisor 
in  the  frequency  plane  of  an  optical  S|)ectrum  analyzer 


Pig.  A2  — Typical  decay  curves  for  light-traiumission  decay  through  a photo- 
dichroic crystal  when  used  in  the  signal  excision  mode 
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Appen'lix  B 

REORIENTATION  CHARACTERISTICS 


To  fully  understand  the  ret-ording  m«H:hanisni  of  the  crystal,  the  reorientation 
characteristics  must  lx-  examined.  Figure*  B1  shows  the  thn'e  pos.sible  orientations  of  the 
defects  and  indicaU-s  the  optical  dipole  moments  of  the  two  absorjition  components 
F^j  and  Fa.^.  Signals  are  recorded  by  inducing  a population  difference  between  orienta- 
tions 1 and  2 in  Fig.  Bl  (populations  Hj  and  ng,  with  the  dimension  cm*2).  Only  these 
two  orientations  are  involved,  since  orientat'on  3 is  symmetric  to  light  incident  from  the 
left. 

To  induce  the  maximum  imbalance  between  Oj  and  02,  the  crystal  must  be  exposed 
to  light  exclusively  in  the  Fa,  or  Faj  bands,  and  the  light  must  be  polarized  parallel  to 
the  appropriato  optical  moment.  First  consider  exposure  to  vertically  polarized  Fa.^  light. 
Defects  in  orientations  2 and  3 will  rotate  into  orientation  1,  and  the  equilibrium  popula- 
tion distribution  will  become 

Hj  = N,  02  = 03  = 0,  (Bl) 

where*  N is  the  total  concentration  of  defects.  The  crystal  is  now  fully  polarizc*d.  (This 
is  an  ideal  case;  in  general  the  maximum  degree  of  alignment  rarely  exceeds  90%. ) If  the 
incident  Fa^  light  polarization  is  now  rotated  45°  so  as  to  have  equal  projections  on  the 
1 and  2 orientations,  the  populations  of  1 and  2 will  equalize,  and  at  equilibrium  the 
polarizing  property  is  lost.  Calculations  with  the  aid  of  Fig.  A1  shows  that  at  equilibrium 
the  {>opulatiun  distribution  becomes 

n,  =n2  = 2N/5,  03  = N/5.  (B2) 

The  distributions  Bl  and  B2  may  be  cycled  indefinitely  by  switching  the  *'’a2 
polarization  between  vertical  and  45°. 


(1)  (2)  (31 


Fig.  Bl  — The  three  pouible  orientations  of  the  defect,  showing  the  corresponding  opti- 
cal dipole  moments  of  the  F^.  and  F;^,  bands.  Absorption  of  a photon  by  any  of  these 
orientations  causes  a rotation  into  one  of  the  other  two  with  an  efficiency  of  0.5. 
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CUinsider  now  tho  same  operation  usinj?  F^i  light,  which  has  only  one  cli|>ole  moment 
jx'r  orientation.  In  this  case  the  population  distribution  b«*comes 

/ij  = 0,  ^2  = ^3  = /V/2.  (B3) 

The  crystal  is  again  fully  aligntnl  (polarized)  but  with  only  half  the  population  difference 
of  the  excitation.  A more  significant  difference  between  the  two  cases  can  b<*  s<*«*n 
when  the  Fa,  light  polarization  is  rotaU‘d  45°.  Examination  of  Fig.  B1  shows  that  at 
equilibrium  the  population  distribution  bt'comes 

«]  = n2  = 0,  n3  = A/.  (B4) 

Although  rij  and  have  been  equilized  and  the  polarizing  property  is  lost,  it  has  been 
done  at  the  exptmse  of  driving  all  defe<’ts  into  orientation  3.  Figun*  B1  shows  that  Fa, 
light  incident  from  the  left  cannot  b«*  absorb«*d  by  orientation  3.  Thus,  to  realign  n, 
rtdative  to  n2  with  Fa,  light,  the  crystal  must  b<‘  first  exposixi  to  Fa^  light.  In  summary, 
ojieration  in  the  F’a.^  band  requires  only  a shift  in  polarization  to  cycle  the  polarizing 
prop«>rty,  but  optiration  in  Fa,  also  requires  an  auxiliary  Fa.^  light  source  to  drive  centers 
out  of  orientation  3.  Of  course  the  latter  could  b<*  any  broilband  light  .source.  These 
population  distributions  are  summerized  in  Table  Bl. 

The  preceding  analysis  deals  only  with  equilibrium  distributions,  but  the  dynamic 
chu^cteristics  can  be  derived,  at  least  for  the  simple  case  of  Fa,  light.  Assume  the 
di.stribution  of  centers  is  random  and  vertically  polarized  light  is  incident  from  the  left 
as  indicated  in  Fig.  Bl.  I.et  I be  the  photon  density  (s'*cm*2),  and  let  A be  the  fraction 
of  these  photons  absorbed  by  the  population  nj  in  orientation  1.  Since  the  quantum 
efficiency  for  reorientation  is  0.5,  the  rate  of  change  of  nj  can  be  written 

n,  =-/A(0.5).  (B5) 

There  is  a probability  of  0.5  that  orientation  2 will  gain  any  defect  lost  by  orientation  1; 
therefore 


J 


/ij  =/A(0.5)(0.5).  (B6) 

The  fraction  of  photons  absorbed  is 

A = l-10''^'.  (B7) 

The  populations  .Tj  and  n2  directly  proportional  to  the  optical  densities  Kj  and  K2, 
which  are  more  easily  observed: 


K,  - fen,, 


(B8) 


where  ft  depimds  on  temperature  and  wavelength  and  has  dimensions  cm^  or  cm^/mJ, 
depending  on  whether  photon  density  or  energy  density  is  measured.  At  514  nm,  ft  is 
given  for  77  K by 


and  for  240  K by 


ft  « 2.00  X 10' cm* 

(B9a) 

ft  - 1.05  X 10'*®  cm2. 

(B9b) 
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Table  B1  — Kquilibrium  Population  Distributions 
of  N Defects  Oriented  as  Shown  in  Fi({.  A1 


Excitation 

Population 

Distribution 

Band 

Polarization 

nilN 

n2/N 

- 

n^/N 

••'a* 

fA, 

fA. 

‘'A, 

Vertical 

45° 

Vertical 

45° 

1 

2/5 

0 

0 

0 

2/5 

1/2 

0 

0 

1/5 

1/2 

0 

Equations  B5  and  B6  now  become 


A:,(E)  = - -|fe/(l-10'^‘)  (BlOa) 

and 

K2(£)  = j 10"'^*).  (BlOb) 

where  E is  the  exposure  symbols.  Equations  BIO  have  the  solutions 

= log  [(10^"  - l)e-''W2  + ij  (Blla) 

and 

^2(0=  |/fo--|K,(f).  (Bllb) 

where  Kq  is  the  initial  optical  density  for  514  nm.  Since  It  is  the  exposure,  Eqs. 

Bll  can  be  written 

^fi(f^)  = logl(10^‘’-l)e-P^  + ll  (B12a) 

and 

K2(E)=^Ko-^Ki(E),  (B12b) 


where  p = k/2  and  has  units  of  cm^/mJ  and  E is  expressed  in  mJ/cm^.  For  514-nm  light 
there  are  2.5  X 10**’  photons/s  per  mW.  Hence  at  77  K,  p = 0.50;  and  at  240  K,  p = 0.25. 


Appendix  C 

TRANSMISSION  CHARACTERISTICS 


C’onsidiT  the  situation  depicted  in  Kig.  Cl.  The  incident  light  on  the  photodichroic 
crystal  consists  of  an  t'  vector  with  components  Oj  and  02  along  the  crystal  axis; 

E = O]  i + 021.  (Cl ) 

If  the  crystal  has  dichroism  (the  absorption  constants  /fj  and  A'2  are  different  along  each 
crystal  axis),  then  the  components  ci|  and  02  will  exp«‘rience  a different  attenuation  as 
they  propagaU*  through  the  crystal,  as  illustrated  in  the  upp<>r  part  of  Fig.  Cl  In  addi- 
tion to  this  effect  there  is  also  a birefringence  induced  by  the  disjx^rsion  effects  of  the 
dichroism.  This  latter  effect  can  be  calculated  from  the  Kronig-Kramers  disfiersion  n*la- 
tion 

2it  r”  (‘*>')  - 

- nn{(jj)  = — I dej  , ((-2) 

* * Cl  * 2 2 

•'0  tu  - w 

where  nj  and  ^2  are  the  indices  of  refraction  along  the  crystal  axis  and  oj  is  the  angular 
frequency.  Computer  calculations  have  shown  that  this  effect  (;an  bt*  apfin-ciable  for 
an  optically  anisotropic  defect  .such  as  the  F^  wnter.  The  dispersion  effe<-ts  art'  additive 
in  the  wavelength  range  between  the  two  absorption  bands  and  subtract  in  the  wings 
The  index -of-rt? fraction  difference  gives  rist*  to  a rotation  of  the  components  O]  and  02 
in  addition  to  their  resp«H;tive  attenuations.  This  phase  angle  is 

2ird(ni  - ) 

6 = , (C3) 

where  d is  the  crystal  thickness.  This  effect  is  illustrated  by  the  lower  part  of  Fig.  Cl. 
Thus  the  E vector  after  propagating  through  the  crystal  can  be  writU-n 

E > '2.1* /2:^  ^-0.434  A2/2-,s/2-,  ^,.4, 


with 


|E(0)|2  » /q  = 2a2  (C5) 

being  the  incident  intensity.  Consider  the  transmission  of  this  E vector  through  the 
analyzer  of  Fig.  Cl.  The  normalized  easy  direction  of  this  analyzer  is 


A 


1 


(i-i). 


(C6) 
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Polarizar  Photodichroic  Analvior 

Crystal 


Pig.  Cl  — Attenuation  and  phase  change  in  a photo- 
dichoric  crystal.  The  upper  part  of  the  diagram 
illustrates  the  relative  attenuation  of  the  amplitude 
components  oi  and  ag  as  they  pass  through  the 
crystal.  Amplitude  sees  an  optical  density  of 
K\I2,  and  ag  sees  /Cg/2.  If  ACg  » AC]  (there  is 
dichroism),  ag  is  greatly  attenuated  compared  to  oi , 
which  effectively  rotates  the  E vector  by  some  angle 
(3  given  by  0 ■ (»r/4)  - 10-<A  i -K  2>/*.  The  lower  part 
of  the  diagram  illustrates  the  phase  change  experi 
enced  by  the  components  O]  and  ag  as  they  pass 
through  the  crystal  due  to  the  birefringence  (given 
by  B2).  In  general  both  effects  are  present 

The  intensity  transmission  of  light  through  the  analyzer  is  then 
/=  |E-A|2 

= sin2  5/2|.  (C7) 


Since  the  intensity  transmission  is  given  by  I/Iq,  this  can  be  written 


